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Abstract. We have developed a spherically symmetric dynamical model of a dwarf 
spheroidal galaxy using the Schwarzschild method. This type of modelling yields con- 
straints both on the total mass distribution (e.g. enclosed mass and scale radius) as well 
as on the orbital structure of the system modelled (e.g. velocity anisotropy). Therefore 
not only can we derive the dark matter content of these systems, but also explore possible 
formation scenarios. Here we present preliminary results for the Sculptor dSph. We find 
that the mass of Sculptor within Ikpc is 8.5 x 10^-*^-^^Mo, its anisotropy profile is tangen- 
tially biased and slightly more isotropic near the center. For an NFW profile, the preferred 
concentration (~15) is compatible with cosmological models. Very cuspy density profiles 
(steeper than NFW) are strongly disfavoured for Sculptor. 



1 Introduction 

The existence of dark matter has been invoked to explain discrepancies in the observed kinematics of 
(systems of) galaxies. Dwarf spheroidal galaxies appear to be one of the most dark matter dominated 
galaxies, with total dynamical mass to stellar light ratios in the order of 100-1000 Mq/L©. The nearby 
dSph galaxies have the additional advantage that individual stars can be resolved, and their red giant 
branch stars are bright enough to measure line-of-sight velocities with errors of a few km s"^ Their 
high dynamical mass-to-light ratios makes these systems ideal to study dark matter halos, especially 
their internal structure and to constrain their inner density profiles. To this end we study the Sculptor 
dwarf spheroidal galaxy using orbit based dynamical models, also known as Schwarzschild models. 



2 Method and data 

In our Schwarzschild models, we assume a specific gravitational potential in which we integrate test 
particle orbits. These orbits are drawn from a dense grid in the space of integrals of motion (in our 
case energy and angular momentum). Using quadratic programming we find the set of weights such 
that the superposition of the weighted orbits matches the data. In our case we fit the light distribution, 
which we assume to be a Plummer sphere, and the velocity dispersion profile. For Sculptor we use the 
2304 line of sight velocities from the datasets by UJ and \2j. In the left panel of Fig. [T] we show the 
velocity dispersion profile we fit. 
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Fig. 1. Left: Velocity dispersion profile for Sculptor. Right: Joint probability density functions for the enclosed 
mass within 1 kpc and the NFW scale parameter (rj for Sculptor. 
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Fig. 2. Left: Resulting anisotropy profiles for Sculptor. Right: Probability density function of the inner slope of 
the dark matter density profile for Sculptor. 

3 Results 

In the right panel of Fig. [T] we show the constraints on the mass and scale radius (r^) in the case an 
assumed NFW |4| profile for Sculptor. The black contours show the 1, 2 and 3cr confidence interval. 
It shows that although M(r<lkpc) is well constrained, the scale parameter (r^) is less so. The concen- 
tration {riml^s ~ 15) implied is compatible with cosmological models. The green contour includes a 
cosmological prior (the mass-concentration relation from |3 |), which gives a slightly better constraint 
on Vs. Since the method finds the distribution function, we can also reconstruct the orbital structure. In 
the left panel of Fig [2] we show the anisotropy profile in the left panel of Fig.[2l The green line shows 
the median value, and the blue regions the 1,2,3 a confidence intervals. The velocity ellipsoid of 
Sculptor is more isotropic in the center compared to the outer regions. We also change the inner slope 
(of) of the DM density distribution p(r) = po(r/r^)"^(l + r/rs)~^ (for NFW, a = I). The probability 
density function of a is shown in the right panel of Fig. |2l Although the slope is not well constrained, 
very cuspy profiles (a > 1.5) seem to be ruled out for Sculptor. More results can be found in ||5l . 
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